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We calculate the interaction between two spherical colloidal particles embedded in the isotropic 
phase of a nematogenic liquid. The surface of the particles induces wetting nematic coronas that 
mediate an elastic interaction. In the weak wetting regime, we obtain exact results for the interaction 
energy and the texture, showing that defects and biaxiality arise, although they are not topologically 
required. We evidence rich behaviors, including the possibility of reversible colloidal aggregation 
and dispersion. Complex anisotropic self-assembled phases might be formed in dense suspensions. 
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Dispersions of small particles or liquid droplets in a 
host fluid, namely colloidal suspensions, are a widespread 
and important state of matter [Q. They are long-lived 
metastable states, of fundamental interest from the point 
of view of collective interactions in complex matter. They 
also have considerable technological importance, e.g., in 
paints, coatings, foods and drugs. To prevent coagula- 
tion due to attractive van der Waals forces, usually the 
particles are treated in order to produce Coulombic or 
steric repulsive interactions. Recently, a novel source of 
repulsion has been reported: the elastic distortion of a 
liquid crystal host Q, i.e., a fluid phase with long-range 
orientational order of the molecules. The repulsion arises 
from the competition between the surface aligning prop- 
erty, which favors, e.g., a radial orientation of the ne- 
matic molecules, and the bulk elasticity, which favors a 
uniform nematic orientation |j| . This effect was shown to 
arise also in other anisotropic fluid hosts, e.g., lyotropic 
solutions of anisotropic micelles B , and in different liquid 
crystal phases, e.g., cholesterics || . Lately, such systems 
were shown to form liquid or solid composite materials 
with unusual properties |@^] . 

It is well known that solid surfaces influence not only 
the orientation of liquid crystals, but also, in general, the 
degree of order of their constituent molecules. In partic- 
ular, a surface can induce a wetting nematic layer even 
above the transition temperature at which the nematic 
becomes an ordinary isotropic liquid [8|-[l0|. Therefore, 
an elastic colloidal stabilization could be achieved also 
above the nematic-isotropic transition. In addition, the 
vicinity of a phase transition may give a critical charac- 
ter to the stabilization mechanism and yield rich phase- 
separation behaviors, as predicted in ]TT[ | for a simpler 
system with a scalar order-parameter. 

To understand these complex collective effects, it is es- 
sential to calculate the interaction between the particles 
precisely. An attempt, based on a quasi-planar approx- 



imation and assuming uniaxiality of the nematic tenso- 
rial order, has been proposed in Ref. Here, we give 
the first exact solution to the problem, including biaxi- 
ality. Our calculation, based on a multipolar expansion, 
is valid for weak surface ordering. We obtain numeri- 
cally the texture between two spherical particles impos- 
ing normal boundary conditions: unexpectedly, we find a 
defect line that is not topologically required and a strong 
biaxiality around it. We give an analytical expression 
for the interaction energy at large distances, which we 
find always attractive. At short distances, comparable 
to the nematic coherence length, we find either attrac- 
tion or repulsion, depending critically on the distance 
to the nematic-isotropic transition. Collective behaviors 
could thus be critically tuned in such systems, in anal- 
ogy with, e.g., the switchable tackiness or wettability of 
liquid crystal polymers [|l3| . 

Nematic liquid crystals are anisotropic liquid phases, 
in which elongated molecules display a long-range ori- 
entational order. Since nematics are non-polar, this or- 
der is described by a symmetric traceless tensorial order- 
parameter Qij (i,j = 1,2,3). The eigenvectors of Qij 
represent the axes of main molecular orientation and its 
eigenvalues describe the amount of ordering and biaxial- 
ity fl4| . Usually, nematics are uniaxial phases, however 
biaxiality naturally arises in inhomogeneous situations, 
e.g., in the vicinity of defects |15|. For a weak induced ne- 
matic order in the isotropic phase, the quadratic Landau- 
de Gennes expansion of the bulk free-energy density E3] 
has the form 



f n ^ Qij Qi : 



1 



-L Qij.kQij.ki 



(1) 



were comma indicates derivation and summation over re- 
peated indices is implied. Here a > (resp. L > 0) 
quantifies the cost of creating (resp. distorting) the ne- 
matic phase. For simplicity, we take a one-constant ap- 
proximation for the gradient terms. We consider spher- 
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ical colloidal particles dispersed in the isotropic phase 
of the nematic. They favor a uniaxial nematic order 
Qij = S$(viVj — \5ij) on their surface, of outward nor- 
mal v>. Here, < So < 1 is the preferred surface scalar or- 
der parameter and Sij the Kronecker delta. At quadratic 
order, the corresponding surface free-energy density is 
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where W measures the anchoring strength. 

Scaling lengths to the nematic correlation length £ = 
(Lja) 1 ! 2 and energies to F = cl£ 3 Sq, the total free- 
energy can be written as 



F = E c 
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where q a = (Q xx - Q yy )/S , qi = Q yz /S , q 2 = Q X z/S a , 

53 = Qxy/So, 94 = Qzz/So, Cq = 1/4, C\ = C 2 = C 3 = 1, 

C4 = 3/4, and w = W/(a£) is a normalized anchoring 
strength. The corresponding equilibrium equations are 
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in the bulk, and 
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on the surface of each colloidal particle. In spherical co- 
ordinates (r p , 9p,4>p) centered on colloidal particle p, the 
general solution of Eq. (|J) that is regular everywhere but 
in r p = can be written as a multipolar expansion: 
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where Yg m (9 p , <p p ) are spherical harmonics and ui{r p ) — 
\/2/(Trr p )Kg +1 / 2 (r p ), in terms of half-integer modified 
Bessel functions. The solution to Eqs. (|J) and (|^) can be 
written as a superposition of a multipolar expansion (|^) 
for each particle, the coefficients of which can be adjusted 
by imposing the boundary conditions (0) for each spher- 
ical harmonic. 

In the following, we consider the case of two identical 
particles, with reduced radius R and centers in z = ±d/2. 
Then, by symmetry, 
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FIG. 1. Order-parameter profile between two spheres of 
reduced radius R = 2 with normalized anchoring strength 
w — 4. The profile has cylindrical symmetry about the axis 
joining the centers of the two spheres. Top: field lines of 
the nematic director n; Middle: contour lines of the uniaxial 
order parameter S; Bottom: contour lines of the biaxial order 
parameter B. The dots indicate the center of the core (the 
biaxial region) of the — | defect Saturn ring. 



where (r, 9, cf>) are spherical coordinates relative to the 
origin, and 
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a(r,9) = ^2^2 aeUe(r p )Pl (cos 9 p ), 
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in which the are the modified Legendre functions 
appearing in the spherical harmonics Yi m - 

Solving for at, fit, je at leading order in the reduced 
separation d, we obtain the asymptotic interaction free- 
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FIG. 2. Eigenvalues of the nematic order parameter Q 
along a diameter of the — i defect Saturn ring whose section 
is shown in Fig. |l]. On the ring's center, Ai and A2 coincide. 

energy (subtracted of the self-energy of each particle) : 

~8irw 2 R 8 exp(-d + 2R) 

" 3d [9 + 3 (3 + w) R + (4 + 3w) R 2 + (1 + 10) R 3 } 2 ' 

(9) 

Hence, the interaction is always asymptotically attrac- 
tive, in disagreement with the estimation of Ref. Jl2| . 

At arbitrary distance d, we solve numerically the linear 
system giving ai, fy, ji by truncating the expansions (||) 
at some order ^ max and checking for convergence. From 
the three eigenvalues A 1 > A2 > A3 (with Ai+A 2 +A 3 = 0) 
we obtain the scalar order-parameter S — §Ai and the 
biaxiality B = ^(A 2 — A3). We call nematic director n 
the direction corresponding to the largest eigenvalue Ai. 

Typical profiles of n, S and B are shown in Fig. [I]. 
The tensor Q is everywhere continuous, however, a de- 
fect (discontinuity) of n appears, in the form of a "Sat- 
urn ring" of strength — ^ . This defect is not topologically 
required: one can expand it to infinity, where it would 
disappear since Q vanishes. Thus, it arises spontaneously 
to minimize the energy, similarly to focal conic defects in 
smectic-A "batonnets" |l(|. Note that our calculation 
does include the defect's core energy, since it resolves the 
gradients of Q also within the core. Our only approxi- 
mation is the neglect of higher-order gradient terms in 
Eq. (]l|): this is however justified as long as £ is much 
larger than the range of molecular interactions. Biax- 
iality develops within the Saturn ring. Again, this is 
surprising, since both the bulk and the surface favor uni- 
axiality. Thus, any modeling of similar systems assuming 
uniaxiality might miss important features. The behavior 
of the eigenvalues of Q along a diameter of the defect 
is shown in Fig. |^. On the ring's center, Ai = A 2 : the 
nematic is locally uniaxial in the direction of the ring, 
however with a discotic-like ordering, since the ring axis 
corresponds to the smallest eigenvalue. 

The behavior of the interaction energy as a function of 
the separation d exhibits many detailed features. At large 



FIG. 3. Diagram of the different types of interaction energy 
profiles. 
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FIG. 4. Dimensionless interaction energy F as a function 
of d for R = 2 and w = 1.5 (continuous line: case 0), w — 2.3 
(dashed line: case la), w = 2.7 (dash-dotted line: case lb). 

d's, it is always attractive, in agreement with the asymp- 
totic expression ([)[). Then, we distinguish six scenarios, 
depending on the value of R and w (Fig. ||). In case 0, 
the interaction is monotonically attractive. In case la 
(resp. 16), a minimum develops, the interaction energy 
at contact being negative (resp. positive), as shown in 
Fig. H In case 2a, 26 and 2c, the interaction energy dis- 
plays first a minimum then a maximum, as d decreases. 
In case 2a, both the maximum and the contact energy 
are negative; in case 26, the maximum is positive while 
the contact energy is negative; finally, in case 2c both the 
maximum and the contact energy are positive (Fig. [5]). 

To understand the collective behavior of such colloids, 
one needs to compare the amplitudes of the barriers and 
of the wells of the interaction energy with k-gT. To sim- 
plify, let us assume that the colloids are treated, e.g., 
with surfactants or polymer coatings, such as to prevent 
aggregation due to van der Waals attraction. Qualita- 
tively, we distinguish three different cases. Whenever the 
interaction energy displays a minimum deeper than k^T, 
attraction occurs, yielding aggregation either at contact 
or at a finite distance. Otherwise, two cases are possible. 
If a maximum is present with a height larger than k^T, 
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FIG. 5. Dimensionless interaction energy F as a function 
of d for 7? = 0.7 and w = 2 (continuous line: case 0), u> = 2.6 
(dashed line: case 2a), w — 3.2 (dash-dotted line: case 26), 
w — 3.6 (dotted line: case 2c). Different scales for F are used 
on the two sides of the vertical line. 
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FIG. 6. Different regimes of the interaction between two 
colloidal particles, as gauged by k^T: I indifferent, R repul- 
sive, A attractive. 1Z is the particle radius in molecular or 
micellar units m, and £ is the nematic coherence length. 



short-range repulsion occurs: the nematic alone could 
stabilize the colloidal dispersion. If, on the other hand, 
the energy is everywhere smaller than k B T, the interac- 
tion is indifferent. Figure || displays the phase behavior of 
colloidal particles of fixed radius when changing the dis- 
tance to the nematic transition. Here we have assumed 
L ~ k B T/m, S$ ~ 0.5, and W ~ k B T/m 2 , where m 
is a molecular or a micellar length. Note that W is a 
microscopic anchoring, stronger than the corresponding 
coarse-grained anchoring that would be measured at a 
macroscopic scale Q . 

In conclusion, we have investigated the interaction be- 
tween colloidal particles dispersed in the isotropic phase 
of a nematogenic liquid. Our calculations show that small 
particles, of size ranging from a few to several hundreds 
nematic correlation lengths, could be reversibly dispersed 
or aggregated by tuning the vicinity to the nematic tran- 
sition, e.g., by changing the temperature. Interesting 
effects might occur when more than two particles are 
brought together. Standard aggregation favors three- 
dimensional close packing. Here, however, the onset of 
extra defects might produce anisotropic self-assembled 
structures although both the embedding fluid and 
the particles are isotropic. Indeed, while it is obvious that 
three particles in a line will create two Saturn rings, it is 
not easy to picture the defects, and the associated energy, 
in the case of three particles forming a triangle, or four 
forming a tetrahedra. Our method allows to calculate 
the textures and free-energies of arbitrary arrangements 
of colloidal particles. Such detailed studies are however 
beyond the scope of this paper, since they require heavy 
numerical simulations: further analyses are required to 
investigate whether gel-like structures (when linear ag- 
gregates are favored) or smectic-like structures (when 
two-dimensional close- packing is favored) could arise [jl9| . 
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